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ABSTRACT 

Local correlation tracking techniques are used to measure proper motions in a series of high angular resolu- 
tion (^0"1) penumbra images. If these motions trace true plasma motions, then we have detected converging 
flows that arrange the plasma in long narrow filaments co-spatial with dark penumbral filaments. Assuming 
that these flows are stationary, the vertical stratification of the atmosphere and the conservation of mass suggest 
downflows in the filaments of the order of 200 m s"' . The association between downflows and dark features may 
be a sign of convection, as it happens with the non-magnetic granulation. Insufficient spatial resolution may 
explain why the estimated vertical velocities are not fast enough to supply the radiative losses of penumbrae. 
Subject headings: convection - Sun: magnetic fields - sunspots 



1. INTRODUCTION 

Sunspot penumbrae appear in the first telescopic ob- 
servations of sunspots made four hundred year s ago (see 
the historical introduction by ICasanova i fT997h . Despite 
this long observational record, we still lack of a phys- 
ically consistent scenario to explain their structure, ori- 
gin and nature. Penumbrae are probably a form of 
convection taking Place in highly inclined strong mag- 
netic fields (Danielson 1961; Thomas & Montesinos "19931 
Igphlichenmaieret al. 1998; Hurlburt et al. 2000; Weiss et aQ 
12004). However, there is no consensus even on this general 
description. For example, the observed vertical velocities do 
not suffic e to transpor t the energy radiated away by penum- 
brae ('e.g..l Spruilll987 . and §|5}, which has been used to argue 
that they are not exclusively a convective phenomenon. The 
difficulties of understanding and modeling penumbrae are al- 
most certainly associated with the small length scale at which 
the relevant physical process take place. This limitation biases 
all observational descriptions, and it also makes the numerical 
modeling challenging and uncertain. 

From an observational point of view, one approaches 
the problem of resolving the physically interesting 
scales by two means. First, assuming the existence 
of unresolved structure when analyzing the data, in 
particular, whe n interpreting t he spectral line asym- 
metries (e.g., Bumba 1960; Grigoriev & Katz 1972; 
Golovkol 11 974: Sanche z Almeida & Lites 1992; Wiehr 
1995t ISol anki & Montavoij Il993t iBellot Rubio 200# 
Sanchez Almeida.i2005a) . Via line-fitting, and with a proper 
modeling, this indirect technique allows us to infer physical 
properties of unresolved structures. On the other hand, one 
gains spatial resolution by directly improving the image 
quality of the observations, which involves both the optical 
quahty of the instrumentation and the application of image 
restoration techniq u es (e.g.. |Mu lie3 fl973a b; Bonet et al. 
Ip^, "2004 , 2005; Stachnik et al. 1983; Lites et al. 1990; 
Title etal. 1993; Sutterlin et al. 2001; Schar mer et al. 2002; 
Rimmelei K004). Eventually, the two approaches have to be 
combined when the relevant length-scales are comparable 
to the photon mean-free-path (see, e.g., |Sanchez Almeida 
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'2001"). The advent of the Swedish Solar Telescope (SST; 
Scharmer et al. 2003a b) has opened up new possibilities 
along the second direct course. Equipped with adaptive 
optics (AO), it allows us to revisit old unsettled issues with 
unprecedented spatial resolution (~0"1), a strategy which of- 
ten brings up new observational results. In this sense the SST 
has akeady discovered a new penumbral structure, na mely, 
dark lanes flanked by two bright filaments ( Scharme r et alJ 
l200llRouPD~ van der Voort et aL2004,) . These dark cores in 
penumbral filaments were neither expected nor theoretically 
predicted, which reflects the gap between our understanding 
and the penumbral phenomenon. 

The purpose of this work is to describe yet another new 
finding arising from SST observations of penumbrae. It turns 
out that the penumbral proper motions diverge away from 
bright filaments and converge toward dark penumbral fila- 
ments. We compute the proper motion velocity field employ- 
ing the local correlation tracking method (LCT) described in 
the next section. Using the mean velocity field computed in 
this way, we follow the evolution of a set of tracers (corks) 
passively advected by the mean velocities. Independently 
of the details of this computation, the corks tend to form 
long narrow filaments that avoid the presence of bright fil- 
aments (Fig. ^p). This is the central finding of the paper, 
whose details uncertainties and consequences are discussed in 
the forthcoming sections. The behavior resembles the flows 
in the non-magnetic Sun associated with t he granulation , 
mesogranulation, and supergranulation (e.g., November 198^ 
^tleetal. 1989; Hirzbergeret al. 1997; Wang et al. 1995| 
Berger et al. 1998). The matter at the base of the photosphere 
moves horizontally from the sources of uprising plasma to the 
sinks in cold downflows. Using this resemblance to gran- 
ular convection, we argue that the observed proper motions 
seem to indicate the existence of downward motions through- 
out penumbrae, and in doing so, they suggest the convective 
nature of the penumbral phenomenon. 

LCT techniques have been ap plied to penumbrae observed 
with lower spatial resolution (see Wang & Zirin 1992; Denke^ 
(1998). Our analysis confirms previous findings of radial mo- 
tions inward or outward depending on the distance to the 
penumbral border In addition, we discover the convergence 
of these radial flows to form long coherent filaments. 

The paper is organized as follows. The observations and 
data analysis are summarized in § |2] The proper motions of 
the small scale penumbral features are discussed in § |3] and 
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§ 13 The vertical velocities to be expected if the proper mo- 
tions trace true mass motions are discussed in § |5l where 
we also consider their potential for convective transport in 
penumbrae. Finally, we elaborate on the implications of our 
finding in § |6l The (in-)dependence of the results on details 
of the algorithm is analyzed in AppendixlAl 

2. OBSERVATIONS AND DATA ANALYSIS 

We employ the original data set of IScharmer et alJ (j^HO^), 
generously offered for p ublic use by the authors. They were 
obtained with the SST ( Scharmer et a n|2Q03ab), a refractor 
with a primary lens of 0.97 m and equipped with AO. The 
data were post-processed to render images near the diffrac- 
tion limit. Specifically, we study the behavior of a penumbra 
in a 28 minutes long sequence with a cadence of 22 s be- 
tween snapshots. The penumbra belongs to the large sunspot 
of the active region NOAA 10030, observed on July 15, 2002, 
close to the solar disk center (16° heliocentric angle). The 
series was processed with Joint Phase-Diverse Speckle (see 
iLofdahl & Scharmer 2003), which provides an angular reso- 
lution of 0."12, close to the diffraction limit of the telescope 
at the working wavelength (G-band, A 4305 A). The field-of- 
view (FOV) is 26" x 40", with pixels 0"041 square. The 
images of the series were corrected for diurnal field rota- 
tion, rigid aligned, destretched, and subsonic Fourier filtered^ 
(modulations larger than 4 km s~' are suppressed). For addi- 
tional details, see Scharmer et al. ( 2002) and the web page set 
up to distribute the data"*. 

Our work is devoted to the penumbra, a region which we 
select by visual inspection of the images. Figure[^ shows the 
full FOV with the penumbra artificially enhanced with respect 
of the umbra and the surrounding sunspot moat. Figure [Q) 
shows the penumbra alone. The qualitative analysis carried 
out in the next sections refers to the lower half of the penum- 
bra, enclosed by a white box in Fig. [^. The results that we 
describe are more pronounced in here, perhaps, because the 
region is not under the influence of a neighbor penumbra out- 
side but close to the upper part of our FOV. The effects of 
considering the full FOV are studied in AppendixlAl 

We compute proper motio ns using the local correlation 
tracking algorithm (LCT) of November & Simon' as 
implemented by Molownv-Horas & Yi ( 1994). It works by 
selecting small sub-images around the same pixel in contigu- 
ous snapshots that are cross-correlated to find the displace- 
ment of best match. The procedure provides a displacement 
or proper motion per time step, which we average in time. 
These mean displacements give the (mean) proper motions 
analyzed in the work. The sub-images are defined using a 2D 
Gaussian window with smooth edges. The size of the window 
must be set according to the size of the structures selected as 
tracers. As a rule of thumb, the size of the window is half the 
size o f the structure that is being tracked (see, e.g., Bonet et al. 
l20Ql . We adopt a window of FWHM 5 pixels (= 0"2), track- 
ing small features of about 10 pixels (= 0"4). The LCT al- 
gorithm restricts the relative displacement between successive 
images to a maximum of 2 pixels. This limit constrains the re- 
liability of the proper motion velocity components, Ux and t/,,, 
to a maximum of ±2 pix per time step, which corresponds to a 
maximum velocity of some 3.8 km s"' . (Here and throughout 
the paper, the symbols Ux and Uy represent the two Cartesian 



components of the proper motions.) Such upper limit fits in 
well the threshold imposed on the original time series, were 
motions larger than 4 km s"' were filtered out by the subsonic 
filter. 

Using the mean velocity field, we track the evolution of 
passively advected tracers (corks) spread out all over the 
penumbra at time equals zero. Thus we construct a cork 
movie (not to be confused with the 28 min long sequence 
of images from which the mean velocity field we inferred, 
which we call time series). The motions are integrated in 
time assuming 22 s per time step, i.e., the cadence of the 
time series. Figure [Q? shows the corks that remain in the 
penumbra after 110 min. (The cork movie can be found in 
[http : / / www . iac . es/proyect /solarhr/penc ork . htmll ) 
Some 30 % of the original corks leave the penumbra toward 
the umbra, the photosphere around the sunspot, or the penum- 
bra outside the FOV. The remaining 70 % are concentrated in 
long narrow filaments which occupy a small fraction of the 
penumbral area, since many corks end up in each single pixel. 
As it happens with the rest of the free parameters used to find 
the filaments, the final time is not critical and it was chosen 
by trial and error as a compromise that yields well defined 
cork filaments. Shorter times lead to fuzzier filaments, since 
the corks do not have enough time to concentrate. Longer 
times erase the filaments because the corks exit the penumbra 
or concentrate in a few sparse points. 

We will compare the position of the cork filaments with 
the position of penumbral filaments in the intensity images. 
Identifying penumbral filaments is not free from ambiguity, 
though. What is regarded as a filament depends on the spatial 
resolution of the observation. (No matter whether the resolu- 
tion is 2" or 0"1, penumbrae do show penumbral filaments. 
Obviously, the filaments appearing with 2" and 0."1 cannot 
correspond to the same structures.) Moreover, being dark or 
bright is a local concept. The bright filaments in a part of 
the penumbra can be darker than the dark filaments elsewhere 
re.g.. iGrossmann-Doerth & Sch midt 1981). Keeping in mind 
these caveats, we use the average intensity along the observed 
time series to define bright and dark filaments, since it has 
the same spatial resolution as the LCT mean velocity map. In 
addition, the local mean intensity of this time average image 
is removed by subtraction of a running box mean of the av- 
erage image. The removal of low spatial frequencies allows 
us to compare bright and dark filaments of different parts of 
the penumbra. The width of the box is set to 41 pixels or 
1"7. (This time averaged and sharpened intensity is the one 
represented in Figs.^) 

The main trends and correlations to be described in the next 
sections are not very sensitive to the actual free parameters 
used to infer them (e.g., those defining the LCT and the lo- 
cal mean intensities). We have had to choose a particular set 
to provide quantitative estimates, but the variations resulting 
from using other sets are examined in AppendixlAl 

A final clarification may be needed. The physical properties 
of the corks forming the filaments will be characterized using 
histograms of physical quantities. We compare histograms at 
the beginning of the cork movie with histograms at the end. 
Except at time equals zero, several corks may coincide in a 
single pixel. In this case the corks are considered indepen- 
dently, so that each pixel contributes to a histogram as many 
times as the number of corks that it contains. 



' The subsonic filter removes fast oscilla tions mostly d ue to p-modes and 
residual jitters stemming from destretching ITitle et alll989B . 
'* |http : / /www ■ solarphysics ■ kva ■ s e/data/l p"02 /| 



3. PROPER MOTIONS 
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We employ the symbol U/, to denote the proper motion vec- 
tor. Its magnitude f//, is given by 

ul = u^+u^, (1) 

with f/j and f/,, the two Cartesian components. Note that these 
components are in a plane perpendicular to the line of sight. 
Since the sunspot is not at the disk center, this plane is not 
exactly parallel to the solar surface. However, the differences 
are insignificant for the kind of qualitative argumentation in 
the paper (see item #0in Appendix|Aj. It is assumed that the 
plane of the proper motions defines the solar surface so that 
the z direction corresponds to the solar radial direction. The 
solid line in Figure |2| shows the histogram of time-averaged 
horizontal velocities considering all the pixels in the penum- 
bra selected for analysis (inside the box in Fig.^). The dotted 
line of the same figure corresponds to the distribution of ve- 
locities for the corks that form the filaments (Fig. ^). The 
typical proper motions in the penumbra are of the order of 
half a km s"'. Specifically, the mean and the standard devia- 
tion of the solid line in Fig.|2]are 0.51 km s"' and 0.42 km s"' , 
respectively. With time, the corks originally spread all over 
the penumbra tend to move toward low horizontal veloci- 
ties. The dotted line in Fig. |2] corresponds to the cork fila- 
ments in Fig. nj), and it is characterized by a mean value of 



0.21 km s"'and a standard deviation of 0.35 km s"'. This 
migration of the histogram toward low velocities is to be ex- 
pected since the large proper motions expel the corks making 
it difficult to form filaments. 

The proper motions are predominantly radial, i.e., parallel 
to the bright and dark penumbral filaments traced by the in- 
tensity. Figure|3]shows the distribution of angles between the 
horizontal gradient of intensity and the velocity. Using the 
symbol / to represent the intensity image, the horizontal gra- 
dient of intensity V/,/ is given in Cartesian coordinates by, 

with the superscript f representing transpose matrix. The in- 
tensity gradients point out the direction perpendicular to the 
intensity filaments. The angle 9 between the local velocity 
and the local gradient of intensity is 

cos 6^=7— — . (3) 

|V,,/| Uh 

The angles computed according to the previous expression 
tend to be around 90° (Fig. |3Jl, meaning that the velocities 
are perpendicular to the intensity gradients and so, parallel 
to the filaments. The radial motions tend to be inward in the 
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Fig. 2. — Histogram of the horizontal velocities in the penumbra as in- 
fen'ed from LCT techniques. The solid line shows the whole set of penum- 
bral points. The dotted line coiresponds to the velocities at the positions of 
the corks after 110 min. The corks drift toward low horizontal velocity re- 
gions. 
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Fig. 3. — Histogram of the angle between the horizontal velocities and 
the horizontal gradients of intensity. They tend to be perpendicular (the his- 
tograms peak at 9 ~ 90°). The same tendency characterizes both all penum- 
bral points (the solid line), and those having cork filaments (the dotted line). 



inner penumbra and outward in the outer penumbra, a system- 
atic behavior that can be inferred from Fig.|3 At the starting 
position of each cork (i.e., the position at time equals zero), 
Fig-Sshows the intensity that the cork reaches by the end cork 
movie (i.e., at time equals 110 min). This forecast intensity 
image has a clear divide in the penumbra. Those points close 
enough to the photosphere around the sunspot will become 
bright, meaning that they exit the penumbra. The rest are dark 
implying that they either remain in the penumbra or move to 
the umbra. Such radial prop er motions are well known in 
the penumbral literature (e.g.. iMulleJllQTSal lDenkedll998t 
[Sobotka et al. 1999). However, on top of this predominantly 
radial flow, there is a small transverse velocity responsible for 
the accumulation of corks in filaments (§E}. 

The corks in Fig.fO' form long chains that avoid bright fil- 
aments and overlie dark filaments. The tracks followed by a 
set of corks finishing up in one of the filaments are plotted in 
Fig. 13 It shows both the gathering at the head of the filament, 
and the tendency to avoid bright features where the narrow 
cork filament is formed. The migration of the corks toward 
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Fig. 4. — Mapping of the final intensities that corks in each location of the 
FOV have by the end of the cork movie. The points in the outer penumbra 
tend to exit the sunspot, and so they are bright in this image. The points in the 
inner penumbra remain inside the sunspot, either within the penumbra or the 
umbra. The white contours outline the boundaries of the penumbra in Fig.fTI 
Points outside the penumbra are also included. 
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Fig. 5 . — Tracks followed by a set of 80 corks that end up in one of the cork 
filaments. The spatial coordinates have the same origin as those in Figs.lTland 
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Fig. 6. — (a) Histograms of the distribution of intensity associated with 
the cork filaments. The original distribution corresponds to corks uniformly 
spread out throughout the penumbra (the solid line). The dotted line repre- 
sents the distribution at the cork filaments. Note the shift. All intensities are 
referred to the local mean intensity, which explains the existence of negative 
abscissae, (b) Same as (a) but in logarithm scale to appreciate the lack of very 
bright features associated with the cork filaments. 



dark penumbral filaments is quantified in Fig. |6p. It shows 
the histogram of intensities associated with the initially uni- 
form distribution of corks throughout the penumbra (the solid 
line), and the final histogram after 110 min (the dotted line). 
A global shift toward dark locations is clear. The change is of 
the order of 20%, as defined by, 

AAi/(T= [M/(110)-/i/(0)]/a/(0)~-0.21, (4) 

with /i/(f) the mean and aj{t) the standard deviation of the 
histogram of intensities at time t min. Figure |6j) contains the 
same histogram as Fig.|6t but in a logarithm scale. It allows 
us to appreciate how the shift of the histogram is particularly 
enhanced in the tail of large intensities. The displacement 
between the two histograms is not larger because the corks do 
not end up in the darkest parts of the dark penumbral filaments 
(see, e.g., Fig.|5}. 

4. FORMATION OF CORK FILAMENTS 

It is important to know which properties of the velocity field 
produce the formation of filaments. Most cork filaments are 
only a few pixels wide (say, from 1 to 3). The filaments are so 
narrow that they seem to trace particular stream lines, i.e., the 
ID path followed by a test particle fed at the starting point of 
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Fig. 7. — The solid line vectors show three velocities of a convergent 
velocity field. The dashed line vector n con'esponds to the vector normal to 
the velocity U;, in the point r, so that ii lJ;,(r) = 0. The dotted line vector 
shows the component of the velocity field parallel to n when moving along 
n. Note that it is anti-parallel to n, which implies that the directional gradient 
Q is negative for convergent field lines. 



the filaments. Then the presence of a filament requires both 
a low value for the velocity in the filament, and a continuous 
source of corks at the starting point. The first property avoids 
the evacuation of the filament during the time span of the cork 
movie, and it is assessed by the results in § |3] Fig. |2l The 
second point allows the flow to collect the many corks that 
trace each filament (a single cork cannot outline a filament.) 
If the cork filaments are formed in this way, then their widths 
are independent of the LCT window width. 

For the corks to gather, the stream lines of different corks 
have to converge. A convergent velocity field has the topology 
shown by the three solid line vectors in Fig.0 As it is repre- 
sented in the figure, the spatial variation of the velocity vector 
U/, in the direction n perpendicular to U/, is anti-parallel to n. 
(U/, n = with |n| = 1.) Consequently, the places where the 
velocities converge are those where f2 < 0, with 



17 = n-[(nV,,)U,,]. 



(5) 



The equation follows from the expression for the variation of 
a vector U;, in t he direction of the vector n, which is given by 
(nV/,)U/, (e.g., ' Bronshtein & SemendvavevlI985L § 4.2.2.8). 
Then the component of this directional derivative in the di- 
rection normal to U/, is given by fl in equation Q- The more 
negative fl the larger the convergence rate of the velocity field. 
(Note that the arbitrary sign of n does not affect 57.) Using 
Cartesian coordinates in a plane, equation (|5} turns out to be. 



Ul dy U} dx 



U,Uy \dU, dU, 
dx dy 



(6) 



The histograms of 51 for all the pixels in the penumbra and for 
the cork filaments are shown in Fig.|8] Convergent and diver- 
gent flows coexist in the penumbra to give a mean Vl close to 
zero (the solid line is characterized by a mean of 4 x 10"^ s"' 
and a standard deviation of 1 .8 x 10"^ s"'). However, the cork 
filaments trace converging flows (the dotted line has a mean 
of-1.1 X 10"^ s"' and a standard deviation of 2.7 x 10"-^ s~'). 
The typical f2 at the corks implies moderate convergent veloc- 
ities, of the order of 100 m s"' for points separated by 100 km. 



5. DISCUSSION 
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Fig. 8. — Histograms of the derivative of tlie velocities in the direction 
perpendicular to the velocity field (Q; see its definition in the main text). As 
in previous figures, the soHd fine corresponds to all the pixels in the penumbra 
whereas the dotted line describes the cork filaments. 



Assume that the observed proper motions trace true mass 
motions. The horizontal velocities should be accompanied by 
vertical motions. In particular, those places traced by the cork 
filaments tend to collect mass that must be transported out of 
the observable layers by vertical motions. The need for mass 
conservation allows us to estimate the magnitude of such ver- 
tical velocities. Mass conservation in a stationary fluid im- 
plies that the divergence of the velocity field times the density 
is zero. This constraint leads to. 



(7) 



dU, ^ dUy 
dx dy 

as proposed by 'November et al." (' 1 987^; 'November 
detailed derivation of the equation is given in Appendix IbI 
The symbol stands for the scale height of the flux of 
mass, which mu st be close to the density scale height (see 
iNovember et al.l 1987; November 1989). We want to stress 
that neither equation (0 assumes the plasma to be incom- 
pressible, nor /ij is the scale height of U,. Actually, the value 
of /i-, including its sign, is mostly set by the vertical stratifica- 
tion of density in the atmosphere (see Appendix B). Figure|9] 
shows histograms of U~ computed using equation ([7) with 
h- =100 km. We adopt this scale height because it is close 
to, but smaller than, the figure measured in the non-magnetic 
Sun by November (1989, 150 km). The density scale hight de- 
creases with temperature, which reduces the penumbral value 
with respect to that in the non-magnetic photosphere. (One 
can readily change to any other h- since it scales all verti- 
cal velocities.) According to Fig.|5] we find no preferred up- 
flows or downflows in the penumbra. The solid line represents 
the histogram of U, considering all penumbral points; it has 
a mean of only 3 x lO""* km s"' with a standard deviation of 
0.39 km s~'. However, the cork filaments prefer downflows. 
The dotted line shows the histogram of U- for the corks at the 
cork filaments. It has a mean of -0.20 km s"'with a standard 
deviation of 0.36 km s~' . 

According to the arguments given above, the cork filaments 
seem to be associated with downflows. The cork filaments 
are also associated with dark features (§|3}. This c ombination 
characterizes the non-magnetic granulation (e.g. Spru it et aP 
[T99O.) . and it reflects the presence of convective motions. The 
question arises as to whether the velocities that we infer can 



transport the radiative flux escaping from penumbrae. Back- 
of-the-envelope estimates yield the following relationship be- 
tween convective energy flux /v, mass density p, specific heat 
at constant pressure cp, and temperature difference between 
upflows and downflows AT, 



(8) 



jSpruillll987l) . Following the arguments by 'Spruit (198'3 
§ 3.5), the penumbral densities and temperature differences 
are similar to those observed in the quiet Sun. Furthermore /v 
is a large fraction of the quiet Sun flux (75 %). Then the verti- 
cal velocities required to account for the penumbral radiative 
losses are of the order of the velo cities in the non-magnetic 
granulation or 1 km s"' (see also iSchlichenmaier & Solankil 
|2003; Sanchez Almeida 2OO51J. One may argue that the 
penumbral vertical velocities inferred above are far too small 
to comply with such needs. However, the spatial resolution 
of our velocity maps is limited. The LCT detects the proper 
motions of structures twice the size of the window, or 0"4 
m our case (see §13. The LCT smears the horizontal veloc- 
ities and by doing so, it smears the vertical velocities too^. 
Could this bias mask large vertical convective motions? We 
believe that it can, as inferred from the following argument. 
When the LCT procedure is applied to normal granulation, it 
leads to vertical velocities of a few hundred m s~', which are 
smaller than the fiducial figure required for equation (|8} to 
account for the quiet Sun radiative flux (1 km s"'). However, 
the quiet Sun radiative flux is deposited in the photosphere by 
convective motions. Therefore a large bias affects the quiet 
Sun estimates of vertical velocities based on LCT, and it is 
reasonable to conjecture that the same bias also affects our 
penumbral estimates. The vertical velocities associated with 
the non-magnetic solar granulation mentioned above have not 
been obtained from the literature. We have not been able to 
find any estimate when the LCT window has a size to track in- 
dividual granules, say, 0."7-0"8^. Then, we carried out an add 
hoc estimate using the quiet region outside the sunspot studied 
by Bonetet al. (2005). The vertical velocities are computed 
employing equation (0 when the FWHM of the LCT window 
equals 0"75 . The inferred vertical velocities ^ have standard 
deviations between 280 m s"' and 110 m s"' depending of the 
time average used to compute the mean velocities (5 min and 
120 min, respectively). One can conclude that the vertical ve- 
locities are some five times smaller that the fiducial 1 km s"' . 
This bias would increase our velocities to values consistent 
with the radiative losses of penumbrae. 

In short, given the limited spatial resolution, the velocities 
inferred from LCT may be underestimating the true vertical 
velocities. If this bias is similar to that affecting the non- 
magnetic granulation, then the observed velocities suffice to 
transport the radiative losses of penumbrae by convection. 

All the discussion above assumes the proper motions to 
trace true motions. However, the proper motion velocities 
disagree with the plasma motions inferred from the Ever- 

^ Equation is linear in Ux and JJy so that it also holds for averages of Ux 
and Uy, rendering averages of I/,; see Appendix IbI 

* November 1 1989) and Hirzberser et al." tl 991^ emp loy a larger window 
to select mesogranulation, whereas Wans et alJ I1995hli do not provide units 
for the divergence of the horizontal velocities. 

^ The finding of low vertical velocities should not depend on the specific 
observation we use. The order of magnitude estimate, equal for all obser- 
vations, leads to 100 m s"' - consider a gradient of horizontal velocities of 
the order of 1 km s"' across a granule 1000 km wide. Equation Q with 
/ij = 100 km renders 0. 1 km s"' . 
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0.12 




Fig. 9. — Histogram of the vertical velocities of the all the pixels in the 
penumbra (the solid line) and only those in the cork filaments (the dotted 
line). Upflows yield positive J/;. The corks seem to be preferentially associ- 
ated with downflows of a few hundred m s"' . 



shed effect, i.e., the intense and predominan tly radial out- 
ward flows deduc ed from Doppler shifts (e. g.. ISolankill2003l 
iThomas & Weis s-'2()Q4). Our proper motions are both inward 
and outward, and only of moderate speed (§|3}. This disagree- 
ment may cast doubts on the vertical velocities computed 
above. Fortunately, the doubts can be cleared up by acknowl- 
edging the existence of horizontal motions not revealed by 
the LCT technique, and then working out the consequences. 
Equation is linear so that different velocity components 
contribute separately to U-. In particular, the Evershed flow 
represents one more component, and it has to be added to 
the proper motion based vertical velocities. The properties 
of the Evershed flow at small spatial scales and large spa- 
tial scales do not modify the velocities in Figure (|9}- On 
the one hand, we apply equation Q to average proper mo- 
tions as inferred with the finite spatial resolution of our ob- 
servations. It is a valid approach which, however, only pro- 
vides average vertical velocities (see Appendix 151. One has 
to consider the contribution of the average Evershed veloci- 
ties, eliminating structures smaller than the spatial resolution 
of the proper motion measurements. On the other hand, the 
large scale structure of the Evershed flow is also inconsequen- 
tial. According to equation 0, a fast but large spatial scale 
flow does not modify U^; add a constant to J/v and f/,., and 
f/j does not change. Only the structure of the Evershed flow 
at intermediate spatial scales needs to be considered, and it 
does not invalidate our conclusion of downflows associated 
with the cork filaments. For the Evershed flow to invalidate 
this association, it would have to provide upflows co-spatial 
with the cork filaments, and so, with dark lanes. However, 
the existence of upflows in dark lanes is not favored by the 
observations of the Evershed effect, which seem to show the 
opposite, i.e., a local correlation between upflows and bright 
lanes (e.g., Beckers & Schroter 1969; Sanchez Almeida et al. 
Il993t ilohan nesson 1993; Schmidt & Schlichenmaier 2000). 
Thus we cannot rule out a bias of the vertical velocities in Fig- 
ure (|9} due to the Evershed flow but, if existing, the observed 
vertical component of the Evershed flow seems to reinforce 
rather than invalidate the relationship between cork filaments 
and downflows. 

6. CONCLUSIONS 



Using local correlation tracking (LCT) techniques, we mea- 
sure mean proper motions in a series of high angular reso- 
lution (0."12) penumbral images obtained with the 1-meter 
Swedish Solar Telescope (SST; Scharmer et al. 2002). Pre- 
vious studies of lower resolution find predominantly radial 
proper motions, a result that we confirm. On top of this trend, 
however, we discover the convergence of the radial flows to 
form long coherent filaments. Motions diverge away from 
bright filaments to converge toward dark filaments. The be- 
havior resembles the flows in the non-magnetic Sun asso- 
ciated with the granulation, where the matter moves hori- 
zontally from the sources of uprising plasma to the sinks in 
cold downflows. Using such similarity, we argue that the 
observed proper motions suggest the existence of downward 
flows throughout the penumbra, and so, they suggest the con- 
vective nature of the penumbral phenomenon. The places 
where the proper motions converge would mark sinks in the 
penumbral convective pattern. 

The presence of this convergent motions is best evidenced 
using tracers passively advected as prescribed by the penum- 
bral proper motion velocity field: see the dots in Fig.^3 and 
Fig. |5] With time, these tracers or corks form filaments that 
avoid the bright features and tend to coincide with dark struc- 
tures. We quantify this tendency by following the time evolu- 
tion of corks originally spread throughout the penumbra. Af- 
ter 110 min, the corks overlie features which are significantly 
fainter than the mean penumbra (the histogram of intensities 
is shifted by 20%; see §|3|l. Assuming that the proper mo- 
tions reflect true stationary plasma motions, the need for mass 
conservation allows us to estimate the vertical velocities at 
the cork filaments, i.e., in those places where the plasma con- 
verges. These vertical velocities tend to be directed downward 
with a mean of the order of 200 m s~' . The estimate is based 
on a number of hypotheses described in detail in ^Sjand Ap- 
pendix O We consider them to be reasonable but the fact 
that the vertical velocities are not direct measurements must 
be borne in mind. The inferred velocities are insufficient for 
the penumbral radiative flux to be transported by convective 
motions, which requires values of the order of 1 km s~' . How- 
ever, the finite spatial resolution leads to underestimating the 
true velocities, a bias whose existence is indicated by various 
results. In particular, the same estimate of vertical velocities 
applied to non-magnetic regions also leads to vertical veloci- 
ties of a few hundred m s"', although we know from numeri- 
cal simulations and observed Doppler shifts that the intrinsic 
granu lar velocities are much larger (e.g.. Stein & Nordlun^ 
1998; Beckers 1981). 

The algorithm used to infer the presence and properties of 
the cork filaments depends on several free parameters, e.g., 
the size of the LCT window, the cadence, the time span of 
the cork movie, and so on. They were originally set by trial 
an error. In order to study the sensitivity of our results on 
them, the computation was repeated many times scanning the 
range of possible free parameters (Appendix lAi. This test 
shows how the presence of converging flows associated with 
dark lanes and downflows is a robust result, which does not 
depend on subtleties of the algorithms. It seems to depend on 
the spatial resolution of the observation, though. 

The downward motions that we find may correspond 
to the ubiquitous downflows indirectly inferred by 
Sanchez Almeida (2005a) from the spectral line asym- 
metries observed in penumbrae. They may also be connected 
with an old observational result by B ecker s & Schrote^ 
il969tl . where they find a local correlation between bright- 
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TABLE A 1 
Summary of uncertainties" 



b 


Case 


Ui, [km s 








1= 


Uz [km s 








start 


end 




start 


end 


start 


end 




reierence^ 


U.jlitU.4z 


ATI _i_n 


-u.z i 




oyit4z 


U.UUltU. jV 


-U.zUltU. JO 




□ 




O ^-7_LO A O 

O.J /±U.4o 


A Tl _LA /I 1 


-0.25 


m _i_/in 


art 1 A ^ 

8y±4z 


U.UUitU.oU 


-U.44±U.DO 




□ 


LCT 0'.'45^ 


0.41±0.31 


0.19di0.20 


-0.18 


91 ±42 


91 ±45 


0.01±0.17 


-0.05±0.08 




§ 


full penumbra' 


0.59±0.50 


0.28±0.38 


-0.16 


90±41 


91 ±43 


A AA 1 A /I O 

U.U0±0.4z 


-0.1d±U.3» 






end time 37 min 




0.26±0.30 


-0.12 




93±45 




-0.14±0.32 




m 


end time 220 min' 




0.11 ±0.23 


-0.28 




83±40 




-0.24±0.26 




m 

§ 


shaipened to I'.' 3^ 
sharpened to 4f'l'' 






-0.20 
-0.20 


91 ±40 
91±41 


90±42 
89±43 










no velocity threshold 


0.58±0.82 


0.22±0.54 


-0.19 


91 ±40 


89±42 


0.00±0.83 


-0.19±0.71 




m 


D = 97 cm' 


0.48±0.38 


0.16±0.23 


-0.11 


90±41 


95±46 


0.00±0.32 


-0.21 ±0.24 




m 


D = 80 cm' 


0.48±0.37 


0.16±0.23 


-0.08 


90±41 


92±46 


0.00±0.31 


-0.22±0.29 






D = 50 cm' 


0.47±0.35 


0.15±0.19 


-0.01 


90±41 


94±47 


0.00±0.30 


-0.24±0.25 






Solar coordinates 


0.52±0.49 


0.26±0.48 








0.01 ±0.40 


-0.20±0.36 






No destretching 


0.55±0.43 


0.20±0.32 


-0.28 


0.91 ±39 


0.91 ±44 


0.00±0.41 


-0.28±0.44 



"The figures after the ± signs do not represent eiTor bars but standard deviations of the histograms. Start and end indicate 
properties at the beginning and the end of the cork movie. The empty entries correspond to values necessarily identical to the 
reference case. 

''Number con'esponding to the list in AppendixlAl 

'^As defined in equation Q. 

''As defined in equation j^. 

"^As defined in equation y). 

'As defined in equation Q. 

S'LCT 0"2 FWHM, penumbra within the box in Fig.^, cork movie ends at 110 min, running box 1"7 wide removed, time 
step 22 s, velocities larger than 3.8 km s"' neglected. 
''FWHM of the window used by the LCT algorithm. 
'Including the portion of penumbra outside the box in Fig.lTk. 
JUsing a time step of 44 s. 

''Width of the running box average subtracted to the intensity image. 

'Diameter of an ideal telescope used to degrade the resolution of the original data. 



ness and Doppler shift with the same sign all over the 
penumbra, and so, corresponding to vertical velocities 
(see also, tS anchez Almeida et al. 1993; Johannesson 1993i 
ISchmidt & Schlichenmaiet. 2000.) . The correlation is similar 
to that characterizing the non-magnetic granulation, which 
also suggests the presence of downflows in penumbrae. Anal- 
yses of spectroscopic and spectropolarimetric sunspot data 
indicate the presence of downflo ws in the external penum- 
bral rim (e.g., Rim mele. .1995: S chhchenmaier & Schmidt! 
p^OO; del Toro Inies ta et all l200lf Bellot Rubio et al. 2003; 
iTYitschler et al.ii2004l) . The association between downflows 
and dark features is particul arly clear in the 0" 2 angular 
resolution SST data studied bv lLanghans et alJj2005l) . Again, 
these downflows may be an spectroscopic counterpart of 
those that we infer However, it should be clear that we also 



find downflows in the inner penumbra, where they do not. 
Whether this fact reflects a true inconsistency or can be easily 
remedied is not yet known. 



Thanks are due to G. Scharmer and the SST team for grant- 
ing access to the data set employed here. G. Scharmer also 
made some valuable comments on the manuscript. The SST 
is operated by the Institute for Solar Physics, Stockholm, at 
the Observatorio del Roque de los Muchachos of the Instituto 
de Astrofisica de Canarias (La Palma, Spain). The work has 
partly been funded by the Spanish Ministry of Science and 
Technology, project AYA2004-05792, as well as by the EC 
contract HPRN-CT-2002-00313. 



APPENDIX 
ROBUSTNESS OF THE RESULTS 

The various free parameters involved in the determination of proper motions were tunned by trial and error to favor the forma- 
tion of narrow chains of corks in the dark penumbral lanes. However, the formation of such chains, their association with dark 
features, and the coincidence with locations of negative horizontal divergence are all robust results. We repeat the computation of 
proper motions and cork filament formation varying the free parameters. The results remain. This section gives a brief account 
of the study. 

We start from the reference case, whose properties have been described in the main text (see also the first row in TableEili. Then 
each one of the parameters characterizing this reference case is varied maintaining the rest unchanged. The variation modifies 
velocities, intensities and angles, a change that we parameterize using the mean and standard deviation of th e corresponding 
histograms. The actual tests are described below. Each one is associated with one or more rows in Table IaTI The number in 
the first column of Table lATI facilitates the identification; it corresponds to the number in the hst below. Such cross-reference is 
needed to follow some of the arguments. 
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1 . We change the size of the LCT window to find out that the smaller the window the clearer the association of cork filaments 
with dark features and downflows. However, such association remains even for windows as large as 0"45 FWHM. 

2. The reference case analyzes the portion of penumbra within the box in Fig.^. When the whole FOV is used, then the cork 
filaments are not so dark and the mean downflows no so intense. However, all trends remain. 

3. We choose the cork filaments as they appear after 1 10 min of evolution of the cork movie. This final time is not critical. 
The qualitative properti es of the reference case are present right from the beginning of the movie, and they are enhanced 
as time increases. Table lA 1 [ includes means and standard deviations at two times, 37 min and 220 min. The latter has been 
computed assuming a time step of 44 s to integrate the cork movie. This fact also allows us to discard any significant effect 
of the time step on the results. 

4. Changing the size of the box used to remove the local mean intensities chan ges t he intensities used in our argumentation. 
However, the dependence is very moderate as attested by the figures in Table lAll 

5. The LCT algorithm gives a few velocities larger than 2 pixels per time step (equivalent to 3.8 km s"'), which is the largest 
displacement used to compute the cross-correlation function (§|2j. These large velocities result from a malfunction of the 
algorithm used to determine the optimum displacement. For this reason the reference case does not include velocities larger 
than 3.8 km s"' . We checked that this threshold does not modify our results in significant way. We repeated the histograms 
for Uh, 0, 1 and U, including velocities larger than 3.8 km s~' . The means of the histograms do not vary by more than 20 % 
(row number|5lin Table lAll . 

6. We investigate the role of the angular resolution of the images used to compute the proper motions and intensities. The 
0."12 angular resolution of the original time series was degraded as if it had been observed with an ideal telescope of 
diameter D. (Each snapshot of the series was convolved with the corresponding Airy disk.) The association between cork 
filaments and dark lanes shows up in the histograms only when the resolution of the degraded images is not far from the 
original one. It is almost gone when D = 50 cm (A^[/a[ ^ -0.007). Table IaTI also includes the cases D = 80 cm and 
D = 97 cm. (Although the SST has D = 97 cm, the original images were restored so that an ideal D = 97 cm telescope 
does indeed reduce the contrast of the small structures used by the LCT algorithm to track.) It seems that the association 
between convergent proper motions and dark lanes depend critically on the angular resolution. This may explain why it 
has remained unnoticed in previous studies based on D = 50 cm class telescopes. 

7. We assume the horizontal and vertical velocities to be horizontal and vertical in a solar coordinate system. However the 
sunspot is slightly out of the disk center so that this assumption is only approximate. In order to evaluate the influence of 
this approximation, the three components of the velocity field were transformed to the true solar horizontal and vertical 
directions. (T wo ro tations are needed.) The transformation does not modify the histograms of velocities in a significant 
way; see Table lAll 

8. The full analysis was repeated using the time series before destretching. The different snapshots were co-aligned assuming 
a rigid shift plus the rotation to be expected from the alt-azimuthal mounting of the SST. This new analysis does not differ 
in a significant way from reference case, so that the histograms seems to be independent of the application of a destretching 
algorithm to the data set. 



DERIVATION OF EQUATION (7) 

This appendix provides a full derivation of equation Q in the vein of the original derivation b v iNo vembeJ il9 89l) . We contribute 
by considering the finite spatial resolution of the observations, and by discussing the validity of the approximation in the context 
of the magnetized penumbral plasma. Such a derivation was found to be necessary by various readers of the original manuscript, 
including the referees. The hypotheses are summarized in the final paragraph. 

The continuity equation for a magnetized plasma with stationary flows is, 

V(pU) = 0, (Bl) 

with the symbols p and U standing for the density and the velocity vector, respectiv ely. The observations provide a kind of 
volume average of the velocities and densities, therefore, in order to apply equation iBl\ to observables, the variables p and 
U have to be replaced with volume averages. We model the volume averages as the convolution with a kernel describing the 
3D region contributing to th e observed signals. Using the property that partial derivatives and convolutions commute (e.g. 
ISanchez Almeida et alJl995 §2), 

(V(pU)) = V(pU)=0, (B2) 

where the angle brackets denote volume average. In general, 

(pU) = (p) (U) + {(p- (p))(U- (U))). (B3) 

We will assume that the cross-correlation between the local fluctuations of density in the resolution element p- (p), and the 
fluctuations of velocity U- (U), are smaller than the product of the mean values, i.e.. 



|((p-(p))(U-(U)))|«|(p)(U)|. 



(B4) 
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Then, 

(pU)~(p)(U), (B5) 

and using equation iB2\ . 

V((p)(U))^0. (B6) 

The inequality JB4t holds if the fluctuations of density are negligible, if the fluctuations of velocity are negligible, or if the 
fluctuations of d ensity and velocity are uncorrected. It is unclear whether any of these hypotheses alone justify the application 
of equation ( IB4t to penumbr ae. H owever, the three of them can cooperate to make the cross correlation negligible, so it is not 

unreasonable to use equation (IB4> as a working hypothesis, even for penumbrae. 

The continuity equation for the actual variables (equation IBID and the mean variables (equation IIB6I ) are formally identical. 
We will use the former for the sake of clarity, dropping the angle brackets from the expressions (U = (U), and p = (p)). However, 
it must be clear that the densities and velocities appearing in all forthcoming equations correspond to the mean velocities and 
mean densities when the actual velocities and densities are averaged over a volume equivalent to the spatial resolution of the 
observations. 

The density of the penumbral plasma presents two distinct types of variation. First, it varies depending on the magnetic field 
and temperature, with overdensities and underdensities that tend to cancel when computing the mean volume averaged density. 
This kind of variation is the only one existing in horizontal planes. Consequently, the horizontal variations of the mean density 
are much smaller than the true horizontal variations of density. Second, there is a systematic drop of density with the height 
in the atmosphere, as required for a plasma trying to be in equilibrium in a strong gravitational field. It is a systematic effect 
affecting all magnetic fields and temperatures and, therefore, it does not cancel in the mean density, which is expected to have a 
significant vertical variation. In mathematical parlance, the mean density at the spatial coordinates x, y and z can be split into two 
components, 

P(x,y,z) = po{z) + pi{x,y,z), (B7) 
with the one describing the systematic vertical variation po larger than the other accounting for the horizontal variations p[. 

Pi < Po- (B8) 
Keeping in mind this condition, the continuity equation can be approximated as, 

V[(po + Pi)U] ~V(poU) = (Vpo)-U + poVU~0. (B9) 



Since pa only vaiies along the vertical direction then. 



and equation ( IB9> can be rewritten as. 



(Vpo)-U=^f/„ (BIO) 

dz 



dU,^dJnp,^ mm^^^ (BID 
oz dz ~ ox ay 

where the symbols Ux, Uy and represent the three Cartesian components of the mean velocity field. With the density scale 
height hp defined as, 

, din On 

/z-;=-— (B12) 



and the vertical velocity scale height given by. 



dz 
ainlf/J 



(B13) 



equation (IBl 1> can be expressed as. 



U,~h:, 



dz ' 

(B14) 



dUx dUy 
dx dy 



where the symbol represents the scale height of the vertical flux of mass, 

/.:' = -^^ = (B15) 

The expression JB14> corresponds to equation Q, which completes its derivation. Note that neither hp nor /!„ have been assumed 
to be constant, therefore, the densities and velocit ies do not necessarily have to vary exponentially. 

An expression formally identical to equation ( IB14> is rea dily d erived from equation ( IBU under the unrealistic assumption 
of constant density. Then VU = 0, which leads to equation (IB14> with h, = h,,. The inference of downflows associated with 
converging flows is based on the assumption h- > (see §|5}- In the unrealistic case of constant density, h^> 0<^ hi, > 0, casting 
serious doubts on the inferred downflows since they follow fron the sign of /i„, which is unknown. However, the density is not 
constant and, consequently, one can have h,>0 with vertical velocities whose magnitude \U~\ decreases with height (/i„ > 0), 
increases with height (/i„ < 0), or is constant (|/i„| oo). Although h^ > is an hypothesis, it seems to be more reasonable 
than the alternative hr < 0. The density is expected to decrease with height {hp > 0) and such a drop is inherited by the flux of 
mass, favoring h^ > 0. In other words, for the assumption h, <0to hold one needs flows of unrealistic amplitude. Consider the 
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following example. According to the arguments in the main text, hp ~ 1 50 km . Assume = -100 km, so that the upflows and 
downflows are reversed with respect to those in Figure|9l Then equation ( IB15> renders hu = -60 km, and so, 

t/,(z) ~ f/,(0) exp(z/60km). (B16) 

This law predicts large unobserved flows as soon as one moves up in the atmosphere. In the mid photosphere with z =150km, an 
upflow of 0.5 km s~' at z = km is amplified to U~(15Q km)=6 km s"^ Such extreme vertical velocities have not been observed 
questioning the underlying assumption h. < 0. Conversely, the positive used in the paper does not cause such problems; 
/ie=100km with hp ~ 150 km lead to h„ = 300km and f/,(150 km)= 0.3 km s-\ 

To sum up, the hypotheses leading to equation Q are: (a) the mean flows are stationary, (b) the cross-correlation between the 
fluctuations of density and velocity in the resolution element are smaller than the product of the mean values, and (c) the mean 
density varies mostly in the vertical direction. 
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